Orientation of Nd"^+ dipoles in yttrium aluminium garnet 
a simple yet accurate model 
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We report an experimental study of the 1.064 /im transition dipoles in neodymium doped yttrium 
aluminium garnet (Nd-YAG) by measuring the coupling constant between two orthogonal modes 
of a laser cavity for different cuts of the YAG gain crystal. We propose a theoretical model in 
which the transition dipoles, slightly elliptic, are oriented along the crystallographic axes. Our 
experimental measurements show a very good quantitative agreement with this model, and predict 
a dipole ellipticity between 2% and 3%. This work provides an experimental evidence for the simple 
description in which transition dipoles and crystallographic axes are coUinear in Nd-YAG (with an 
accuracy better than 1 deg), a point that has been discussed for years. 

PACS numbers: 42.70.Hj, 42.55.Rz 



While Nd-YAG is one of the most (if not the most) 
commonly used solid-state laser crystals, the exact ori- 
entation of transition dipoles within it is still, paradoxi- 
cally, an unresolved problem. This is probably owing to 
the fact that for most applications it is sufficient to con- 
sider the Nd-YAG crystal (usually grown along the (111) 
crystallographic axis) as isotropic, although it has been 
known for long [11,13 thatNd3+ ions in this configuration 
rather see a D2 symmetry, with six possible dodecahedral 
orientations. The influence of crystal symmetry on dipole 
orientations has been previously studied in saturable ab- 
sorbers such as Cr-YAG and Tm-YAG 0,0. In 
the first case, it has been clearly established that transi- 
tion dipoles were aligned with the crystal axes (labeled 
(100), (010) and (001)) [li] while in the second case it 
has been found that they were rather coUinear with the 
(110), (Oil) and (101) directions 0. In the case of Nd- 
YAG, the answer to this question is still unclear in spite 
of several previous studies involving in particular dynam- 
ical polarization effects in Nd-YAG lasers [1, [l3, 11 1 . 



In this paper, we propose a new approach to probe the 
orientation of transition dipoles in Nd-YAG, by measur- 
ing the coupling constant between two linearly-polarized 
orthogonal modes of a laser cavity for different cuts of the 
gain crystal, using a steady-state method similar to the 
one described in |12] . The measured coupling constant is 
a dimensionless ratio between cross-saturation and self- 
saturation coefficients, which is relatively independent of 
most laser parameters (pumping rate, birefringence,. . .), 
hence a good indicator for testing the validity of theoret- 
ical models. Our study deals for the most part with the 
1064.15 nm emission line, sometimes referred to as R2, 
between the upper doublet of "^^3/2 and the Y3 level of 
'*/ii/2. As a matter of fact, it is known from previous 



at low pumping rates. 

The paper is organized as follows. We first propose 
a theoretical model for calculating the coupling constant 
between two orthogonal modes of a Nd-YAG laser cavity, 
on the assumption that transition dipoles are coUinear 
with YAG crystallographic axes. Starting from the very 
simple case of linear transition dipoles, we then generalize 
the model to the case of (slightly) elliptical dipoles, show- 
ing a very good agreement with the experimental value 
of the coupling constant published in 'l2!| with (lll)-cut 
Nd-YAG. We then apply this model to the description of 
our own experimental configuration, namely a two-mode 
laser cavity using as a gain medium a Nd-YAG crystal 
either (111) or (lOO)-cut. In both cases, the measured 
value of the coupling constant is compared with the cor- 
responding theoretical prediction. The results, implica- 
tions and perspectives of this work are finally discussed. 



studies [13|, |14| that the Rl line at 1064.4 nm (between 



the lower doublet of ^-^3/2 and the Y2 level of '*/ii/2) has 
a very small contribution to the overall gain, especially 



EXPRESSION OF THE COUPLING CONSTANT 
IN THE CASE OF LINEAR DIPOLES 



In the usual semiclassical description of lasers, the 
transition dipole is often modeled by an operator d as- 
sociated with the linear vector d = du, whose interac- 
tion with a linearly-polarized electric field E = Ex is 
described quantum mechanically by the dipolar Hamilto- 
nian — d-E. Assuming furthermore that dipole coherence 
lifetime is much shorter than photon lifetime r and pop- 
ulation inversion lifetime Ti (which is indeed the case in 
Nd-YAG lasers [l3|)i the dipolar interaction can be de- 
scribed by the following terms in the rate equations for 
the population inversion density N and field intensity / 
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= -T^T 1-COS (X,U) 
J-1 

acNI cos^ (xTu) , 



(1) 



where Is is the saturation intensity, a the interaction 
cross section, c the speed of Ught in vacuum, and where 
u and X are assumed to be unit vectors. The overall rate 
equations for / and in a simple Lamb's laser model 
then read Il5l : 
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where W is the pumping rate. As regards Nd-YAG lasers, 
the latter equations are sufficient to describe satisfacto- 
rily most of the experimentally observed phenomena such 
as relaxation oscillations [lij or spiking during laser turn- 
on [13]. However, they do not allow an accurate descrip- 
tion of mode coupling in a Nd-YAG laser cavity [l^ . 



General expression of the coupling constant 



regime and for near-threshold operation (i.e. Ij/P^ ^ 1), 
equations ([3]) and ([4]) can be rewritten as : 
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and 
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(6) 



i=l 



The coupling constant C, initially defined by Lamb as 
the ratio between cross-saturation coefficients and self- 
saturation coefficients ^is"] in the case of lasers with short 
population inversion lifetime (Ti ^ r), can be general- 
ized to other kinds of lasers like Nd-YAG by the following 
(more general) definition, involving small variations from 
the steady-state regime [12] : 
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(A/i/A72)(A/2/A7i) 
(A/i/A7i)(A/2/A72) 



(7) 



Using this definition, a straightforward calculation start- 
ing from equations ([6]) and using equations ([5|) for the 
expression of Ni leads to the following formula : 



To this end, one must take into account, in addition 
with the two laser modes, the existence of several possible 
orientations for the transition dipoles. In the following 
model, we shall assume three possible orientations corre- 
sponding to the unitary vectors Ui, U2 and U3, and as- 
sociated with population inversion densities Ni,Nn and 
N^. We furthermore consider, in keeping with |12l ] and 
with the experiment described later on in this paper, the 
case of a laser with two modes linearly polarized along 
the xi and X2 unitary vectors, and associated with in- 
tensities Ii and l2- The semiclassical equations for the 
dipolar interaction ^ can be generalized as follows : 
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C = 
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(8) 



One important point is that this expression depends only 
on the overall geometry, making the coupling constant a 
useful tool for studying dipoles orientation. 



Application to the case of the (111) crystal 
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where we have introduced two possibly different satura- 
tion intensities and . The overall rate equations ([2|) 
for Ni and Ij become in this case : 



diV, 



diV,: 



N, 

W -- 

dt Ti ' dt 



and —r- 

dt 



I. 



dl^ 
dt 



, (4) 



int 



where we have introduced different loss coefficients for 
each mode = 1/tj. In equations W has been 
chosen to be independent of i, which corresponds to the 
case of isotropic pumping (this has been checked experi- 
mentally, see further in this paper). In the steady-state 



The case of the (111) crystal is by far the most common 
for Nd-YAG lasers. In this configuration, the laser wave- 
vector k is along the (111) axis, while the laser electric 
field lies in the transverse plane. Let us assume that 
the two laser modes are linearly polarized along the two 
following transverse unitary vectors : 



Xl 




and X2 = — 7= 

V6 



(9) 



where the coordinates are expressed in the base of the 
crystallographic axes. Although this orthogonal base of 
the transverse plane has been arbitrarily chosen, it should 
be pointed out that this choice does not affect the final 
expression for the coupling constant (calculation shown 
in the appendix). Differently speaking, the coupling con- 
stant is independent, in the (111) case, of the crystal 
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Figure 1: Comparison between usual linear dipole cou- 
pling (left) and elliptical dipole coupling (right). The latter 
has been phenomenologically introduced to allow for cross- 
coupling between dipoles from different sites. 



orientation. Following the authors of 0, we consider 
transition dipoles dui, du2 and du^ coUinear with the 
crystal axes, namely : 
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and U3 
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This hypothesis will be self-consistently confirmed by 
our experimental results later on in this paper. Expres- 
sion ([8]) immediately leads in this case to the coupling 
constant value C = 1/9 ~ 0.11. 

This value can be compared with the experimental 
measurement C ~ 0.16 ± 0.03 from reference fl5|. The 
discrepancy is attributed to the fact that dipoles from 
each site are not perfectly linear, but rather slightly 
elliptic, as will be described in what follows. 



Case of the (111) crystal 

The latter equations can be applied to the previously- 
studied case of the (111) Nd-YAG crystal, using the vec- 
tors defined in ([9]) and (fTO|) . This leads to the following 
rate equations for the population inversion densities : 
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and for the laser modes's intensities : 
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We eventually obtain the following expression for the cou- 
pling constant in the presence of small elliptical dipolar 
coupling, up to the first order in /3 : 



(11) 



Based on the experimental result C ~ 0.16 ± 0.03 of 
reference [13], the estimate /3 ~ 2.75% (±0.5%) can be 
deduced from expression (jlip . A possible explanation 
for this ellipticity id the existence of arbitrarily-oriented 
Nd'^+ ions residing in defect sites of the YAG matrix, 
that would induce energy transfer between dipoles from 
different intrinsic crystal sites. 



TAKING INTO ACCOUNT CROSS-COUPLINGS 
BETWEEN DIPOLES FROM DIFFERENT 
CRYSTAL SITES 

In order to account for possible cross-couplings be- 
tween dipoles from different crystal sites, we shall assume 
a small dipole ellipticity /3 <C 1. In this phenomcnolog- 
ical description, illustrated on figure [U the field-dipole 
interaction ([3]) takes the new following form : 
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= aclj ^ [cos^ (xTTui) + P sin^ (xTTui)] • 



This model will be used to calculate a new expression for 
the coupling constant in both the (111) and (100) cases. 



Case of the (100) crystal 

We now turn to the case of a (100) crystal, on the rea- 
sonable assumption that dipole ellipticity /? is indepen- 
dent of crystal cut. In this new configuration, the laser 
wave-vector k is aligned with one crystallographic axis 
(say U3) and the two orthogonal laser modes are linearly 
polarized along two transverse axes defined by : 



and X2 = 



where the coordinates have been expressed in the base 
(ui,U2,U3) of the crystallographic axes. With the lat- 
ter definition, the polarization directions of the two laser 
modes xi and X2 make an angle a with the crystal (or 
dipoles) axes Ui and U2 . This leads, up to the first order 





in f3, to the following rate equations for the p 
inversion densities : 
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where the following notations have been introc 

A = cos^ a+/3sin^ a and B — sin^ a+Pcof 

Similarly, one gets the following rate equatio: 
modes's intensities : 



— ± =ac [NiA + N2B] h , 




— ^ = (Tc[NiB + N2A]l2 . 



It should be mentioned that the influence c 
been neglected in this analysis since it would in 
terms on the order of 0^ or smaller. One final 
C = AA^B^/iA^ + ^2)2, which reduces, up b 
order in and /3, to the following expression 

This expression will be used later on in this 
compare this theoretical model with data 
experiment. 



EXPERIMENT 

We now turn to the description of our exj 
setup, sketched in figured! We use a 18-cm-L 
laser cavity with a 2.5-cm-long Nd-YAG crj 
gain medium. The cavity also contains a IC 
uniaxial birefringent crystal (YVO4) cut at 
its optical axis, in order to spatially separat 
orthogonal modes of the cavity. We have cht 
the two perpendicular ordinary and extraordin 
izations correspond to the two spatially separated modes 
inside the cavity and that no significant cross-coupling 
inducing forked eigenstate operation [l9| was induced 
inside the Nd-YAG crystal (which imposes a specific 
crystal orientation in the (100) case, see further). We 
have also checked, using a Fabry-Perot analyzer, that 
each one of the two perpendicular laser modes was 
longitudinally single-mode. This is probably due to the 
fact that the YVO4 crystal acts as an etalon and also 
due to the near-threshold operation of the laser. 

The razor blades, placed in the vicinity of the two 
separated beams, are intended for creating additional 
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Figure 2: Sketch of our experimental leiser cavity setup. Each 
razor blade can be translatory moved perpendicularly to the 
laser axis, in order to create additional losses to the corre- 
sponding laser mode. Ml and M2 are the cavity mirrors. 



losses to the corresponding laser modes (i.e. changing 
7i and 72). We limit the losses introduced by the razor 
blade to about 1%, allowing us to neglect the beam 
truncation and to keep a good overlap of the beams in 
the active medium. The intensities of the two orthogonal 
modes are then monitored on two photodiodes after 
being separated by a polarizing beam splitter. Typically, 
the position of one razor blade is changed periodically in 
time (using a piezoelectric transducer), and the relative 
dependence of Ii and I2 is monitored on an oscilloscope. 
We successively introduce razor blade losses to the 1 
and 2 modes, thus obtaining experimental values for 
(A/2/A7i)/(A/i/A7i) and (A/i/A72)/(A/2/A72). 
The coupling constant is eventually deduced from 
expression ([7]). 

It is worth noticing that this result is independent, at 
least in the framework of our theoretical model, of most 
laser parameters, in particular saturation intensities 
and /g, pumping rate and calibration of photodiodes. 
The dominant error source in this measurement is inten- 
sity self-modulation, which introduces an uncertainty in 
the slope measurement. Another error source is resid- 
ual pumping anisotropy. In order to assess the contribu- 
tion of the latter effect, we have checked experimentally 
that pumping light had no preferential polarization axis 
with an accuracy better than a few percents, resulting in 
a relative error of a few percents on the measurement 
of C. The overall measurement error is estimated to 
SC/C -10%. 



Case of the (111) Nd-YAG 

In a first experiment, we use a typical (111) Nd-YAG 
crystal as a gain medium, with arbitrary orientation. The 
raw output light from the fibred laser diode is focused on 
the crystal for optical pumping. Following the method 
described above for the measurement of the coupling con- 
stant, we have obtained the data reported on figure [3l 
This leads to the following value : C~0.15±0.015, in 
very good agreement with the previously-published value 
C ~ 0.16±0.03 ^l3|. This agreement is probably owing to 
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Figure 3: Experimental curves showing I2 versus /i in the 
(111) case when the loss rates 71 (left curve) and 72 (right 
curve) are changed. Slopes measurements provide the fol- 
lowing values : (A/2/A7i)/(A7i/A7i) = -0.3 ± 0.03 (left 
curve) and (Ah / AJ2) / {Ah / A'y2) = -0.5 ± 0.05, leading to 
n = 0.15 + 0.ni5. 
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Figure 4: Relative orientation, around the longitudinal cavity 
axis, of the YVO4 crystal (left) and the (100) Nd-YAG crys- 
tal (right). As can be seen on this sketch, the YAG crystal 
has been oriented in order to align its residual birefringence 
axes with the YVO4 crystal axes, in order to ensure that the 
cavity eigenmodes are linearly polarized along the YVO4 axes 
(instead of being forked modes) . 

the fact that the coupling constant is relatively indepen- 
dent of most laser parameters, as pointed out previously, 
provided the two cavity modes are properly separated by 
the YVO4 crystal and spatial hole burning can be ne- 
glected (which can be shown to be the case here and in 
the work of reference [l^ ) . Taking into account this new 
experimental value (with smaller error bars than in [l^). 
the following finer estimate of the ellipticity f3 can be 
deduced from equation ([IT]) : /? ~ 2.2% (±0.2%). 



Case of the (100) crystal 

In a second experiment, the (111) gain medium is 
replaced by a (100) Nd-YAG crystal, manufactured by 
the German company FEE GmbH. The crystallographic 
axes are known precisely by X-ray analysis. A small 
birefringence is observed for this crystal, probably due 
to mechanical stress from the mount. In order to avoid 
the appearance of forked modes in the laser cavity [l^ , 
it is necessary to properly align the birefringence axis of 
the YAG crystal with the ordinary and extraordinary 
polarization axis of the YVO4 crystal. As illustrated 
on figure 31 such a configuration is experimentally 



Figure 5: Experimental curves showing h versus h in 
the (100) case (angle between laser and crystal axes ; 
10±1 deg)when the loss rates 71 (left curve) and 72 (right 
curve) are changed. Slopes measurements provide the fol- 
lowing values : (A72/A7i)/(A/i/A7i) = -0.11 ±0.01 (left 
curve) and (A/i/A72)/(A/2/A72) = -0.1 ± 0.01, leading to 
C = 0.011 ±0.001. 

obtained, on our setup, with an angle a = 10 ± 1 deg 
between the crystal axes and the YVO4 axes. With 
such an alignment, the cavity eigenmodes are lin- 
early polarized and coincide with the ordinary and 
extraordinary polarization axis of the YVO4 crystal, 
making the experiment suitable for comparison with 
our theoretical model. It is then possible to measure 
the coupling constant between both orthogonal modes 
following the previously-described method, as reported 
on figure [H The result is C ~ 0.011 ± 0.001, to 
be compared with the theoretical value predicted by 
our model (equation (jl3p with a = 10 ± 1 deg and 
13 ~ 2.2% ± 0.2%), namely C ~ 0.011 ± 0.001. This 
remarkable agreement is an experimental evidence for 
the simple theoretical description proposed in this paper. 



CONCLUSION 

To summarize, we have proposed a simple yet ac- 
curate theoretical model for describing the orientation 
of transition dipoles in a Nd-YAG crystal, by making 
the following hypotheses : first, dipoles are slightly 
elliptic; second, their main directions are coUinear with 
the crystallographic axes. This model has been tested 
by making experimental measurements of the coupling 
constant between two orthogonal modes of a Nd-YAG 
laser for two different crystal cuts ((111) and (100)). A 
remarkable quantitative agreement between theory and 
experiment has been observed. In particular, this study 
is an experimental evidence for the fact that transition 
dipoles are indeed aligned with crystallographic axes of 
the Nd-YAG crystal, with an accuracy better than 1 deg. 

The significant reduction of coupling between orthogo- 
nally polarized modes in a laser cavity using a (lOO)-cut 
gain medium, as demonstrated in this paper, could be 
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used to significantly increase the stability of bi-frequency 
lasers, with possible applications in the field of lidar or 
multioscillator ring laser gyroscopes. Our work predicts 
that the more favorable situation for such applications 
will occur when the crystallographic axes are aligned 
with the laser cavity axes, with a minimum achievable 
coupling constant as small as 4/3^ ~ 2.5 10"'^. Further- 
more, the simple and original protocol proposed in this 
paper could be applied to other kinds of solid-state gain 
media or saturable absorbers, in order to probe the 
orientation of their active dipoles. 

The authors are happy to thank Philippe Goldner from 
Chimie ParisTech and Daniel Rytz from FEE GmbH for 
helpful discussions. 
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It is a remarkable fact that expressions ([16 




are 



independent of a. Using expression ([8]) for the coupling 
constant eventually leads to C — 1/9, independently of 
the orientation between the (111) crystal and the laser 
modes. 



APPENDIX: CALCULATION OF THE 
COUPLING CONSTANT FOR AN ARBITRARY 
CRYSTAL ORIENTATION IN THE (111) CASE 
WITHOUT DIPOLE ELLIPTICITY 

In this appendix, we will show that in the case of a 
(111) crystal with linear dipoles oriented along the crys- 
tallographic axes, the coupling constant is independent of 
the orientation between the crystal axes and the direc- 
tions of the laser modes. To this end, we consider linear 
dipoles oriented along the crystallographic axes, namely : 



(14) 



Any pair of orthogonal directions for the laser modes in 
the transverse plane can be obtained by rotating (by the 
appropriate angle a) the initial base defined by ([9|) , which 
reads : 
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A straightforward although tedious calculation leads, 
from equations (|14p and (|15p . to the following expres- 
sion : 
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